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I.  INTRODUCTION 


A  series  of  nonperforating  ballistic  impact  experiments^  were  conducted  by  the  US  Army 
Combat  Systems  Testing  Agency  (USACSTA)  to  characterize  the  ballistic  shock  environment  in 
rolled  homogeneous  armor  (rha)  targets.  One  of  these  experiments  involved  the  normal  impact 
of  a  20  mm  fragment  simulator  projectile  (fsp)  on  914  mm  x  914  mm  x  70  mm  rha  plate  at  an 
impact  velocity  of  1508  m/s.  Normal  displacements  were  measured  at  radial  distances  of  105 
mm  and  241  mm  from  the  center  of  the  plates,  and  radial  strains  were  measured  at  radial 
distances  of  100  mm,  170  mm  and  240  mm  from  the  center  of  the  plates.  The  EPIC>2  hydro- 
code^  was  used  by  the  US  Army  Ballistic  Research  Laboratory  to  model  these  experiments  and  to 
calculate  the  plate’s  responses  to  impact  of  the  fsp.  Comparison^  of  the  measured  and  EPIC-2 
calculated  responses  (normal  displacement  and  radial  strain)  showed  that  the  EPIC-2  could 
calculate  the  response  of  projectile-impacted  targets  within  the  limitations  of  the  code. 

EPIC-2  is  an  axisymmetric  code;  therefore,  it  is  limited  to  axisymmetric  targets  and  to 
normal  impacting  projectiles.  If,  however,  the  EPIC-2  code  is  modified  to  calculate  the  projectile 
loading  history  of  the  target,  this  loading  history  can  be  used  as  input  for  a  finite  element 
analysis  of  ballistic  shock  in  normal  impacted,  non- axisymmetric  targets.  This  report  describes 
the  EPIC-2  calculated  loading  history  of  a  20  mm  fsp  with  an  impact  velocity  of  1508  m/s  on  a 
70  mm  thick  rha  plate,  the  application  of  this  loading  history  in  an  ADINA^  finite  element 
analysis  of  the  plate,  and  the  rei^'ilts  of  the  comparison  of  the  EPIC-2  and  ADINA  calculations. 


n.  EPIC-2  CODE  AND  CALCULATED  LOADING  fflSTORY 


The  major  characteristics  of  the  EPIC-2  code  are  listed  in  Table  1.  The  code  performs 
elastic- plastic  impact  computations  in  two  dimensions  for  axisymmetric  and  plane  strain 
problems  with  either  free  or  fixed  boundaries.  It  is  based  on  a  Lagrangian  finite  element  and 
lumped  mass  formulation  in  which  the  equations  of  motion  are  integrated  directly.  Nonlinear 
material  strength  and  compressibility  efi'ects  are  included  to  account  for  elastic- plastic  flow  and 
wave  propagation.  The  code  has  material  descriptions  for  strain  hardening,  strain  rate  effects, 
thermal  softening  and  failure.  It  uses  a  constant  strain,  triangular  finite  element  which  is  well 
suited  to  represent  the  severe  distortions  occurring  during  high  velocity  impact. 

The  use  of  an  axisymmetric  code  to  model  the  USACSTA  experiments  was  justified  since 
the  impact  point  in  the  experiments  was  the  center  of  the  plate.  EPIC-2  calculations^  had  shown 
that  the  plate’s  response  would  be  axisymmetric  and  independent  of  the  plate’s  geometry  and 
edge  boundary  conditions  for  a  period  of  time  after  impact. 

One  of  the  variables  calculated  by  the  EPIC-2  code  for  each  computational  tin;e  interval 
was  the  total  linear  momentum  of  the  plate.  The  code  was  modified  to  calculate  the  loading 
history  of  the  projectile  on  the  plate  by  calculating  the  time  rate-of-change  in  the  plate’s  total 
linear  momentum  over  each  computational  time  interval.  Although  the  computational  time 
interval  for  the  impact  velocity  was  in  the  order  of  0.05  fts,  the  loading  was  printed  out  at  1.0  ^ 
intervals. 
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Table  1.  EPIC-2  Code  Characteristics 

1 _ 

DISCRETIZATION: 

FINITE  ELEMENT  METHOD 

•  2D  constant  strain  triangles 

•  Lumped  mass  formulation 

MESH  DESCRIPTION: 

LAGRANGIAN 

MATERIAL  MODEL: 

CONSTITUTIVE  MODEL 

•  Incremental  elastic-plastic 

•  Von  Mises  yield  criterion 

•  Compressibility  efl’ects 

•  Strain  rate  effects 

•  Strain  hardening 

•  Thermal  softening 

EQUATION  OF  STATE 

•  Mie-Gruneisen 

FAILURE  CRITERIA: 

VOLUMETRIC  STRAIN 

EFFECTIVE  PLASTIC  STRAIN 

POST-FAILURE  MODELS: 

PRESSURE  CUTOFF 

SHEAR  AND  TENSION  FAJLURE 
TOTAL  FAILURE 

Figure  1  shows  the  computed  loading  history  of  the  plate  for  the  impact  velocity  of  1508 
m/s.  The  loading  has  a  duration  time  of  68  ^^a  and  a  maximum  peak  value  of  6.3  MN.  For  the 
first  30  fis  the  loading  is  very  oscillatory  with  differences  of  up  to  5  MN  between  maximum  and 
minimum  load  values.  There  exists  at  46  /is  a  negative  loading,  implying  that  the  projectile  is 
pulling  the  plate  rather  pushing  against  it.  A  more  reasonable  explanation  is  in  terms  of  the 
dilatation  wave  produced  by  the  impacting  projectile.  The  time  at  which  the  negative  loading 
occurs  is  the  arrival  time  of  the  wave  at  the  impact  point  after  the  wave’s  second  round  trip 
through  the  70  mm  thickness  of  the  plate.  Since  this  is  a  tension  wave  when  it  arrives,  the  plate 
experiences  a  pulling  rather  than  a  pushing. 
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Force  (MN) 


m.  ADINA  CALCULATIONS  AND  EPIC-2  COMPARISONS 


The  ADINA  finite  element  model  of  the  plate  is  shown  in  Figure  2  and  it  is  identical  to 
the  EPIC-2  model.  The  plate  is  modeled  as  a  circular  one,  914  mm  in  diameter  and  with  free 
boundaries,  and  it  consists  of  1029  nodes  and  1760  triangular  elements.  The  horizontal  and  the 
vertical  distances  between  the  nodes  are  7.62  mm  and  8.67  mm,  respectively.  Only  the  first  160 
mm  of  the  plate’s  457  mm  radius  is  shown.  One  third  of  the  loading  shown  in  Figure  1  is  applied 
to  each  of  the  top  first  three  radial  nodes  as  indicated  in  Figure  2.  (The  20  mm  fsp  used  in  this 
experiment  was  completely  destroyed.  However,  recovery  of  similar  fsp’s  after  lower  impact 
velocity  experiments  showed  the  contact  surfaice  of  the  projectile  mushroomed  from  a  20  mm 
diameter  to  a  33  mm  diameter.) 

The  analysis  is  nonlinear,  small  deflection  and  small  strain.  The  time  integration  scheme 
used  is  the  Newmark  method  with  a  time  step  increment  of  1  /xs,  and  the  equilibrium  iteration 
method  used  is  the  full  Newton  with  line  search.  The  material  model  of  the  plate  is  bilinear 
elastic-plastic  in  which  Young’s  modulus  E  =  145.1  GPa,  Poisson’s  ratio  u  =  0.275,  yield  stress 
(Ty  =  916.2  MPa,  and  strain  hardening  modulus  Ej.  =  510.1  MPa;  and  the  density  of  the  plate  is 
7876  kg/m®  The  same  material  properties  are  used  in  the  EPIC-2  calculations. 
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e  2.  ADINA  finite  element  model  of  70  mm  thick  plate. 


Figure  3  through  9  are  plots  of  both  the  ADINA  and  EPIC-2  calculated  responses  at  the 
same  locations  on  the  back  of  the  plate.  Figure  3  shows  excellent  agreement  of  the  normal 
displacements  at  106  mm,  and  Figure  4  shows  good  agreement  of  the  normal  displacements  at 
244  mm.  For  time  t  >  lOO  /js,  the  ADINA  displacements  are  less  than  the  EPIC-2  displacements 
at  244  mm.  Figure  5  shows  good  agreement  of  the  radial  strains  at  100  mm  except  for  50  /zs  <  t 
<  75  fia.  During  this  time  interval  the  EPIC-2  maximum  positive  strain  is  approximately  three 
times  greater  than  the  ADINA  maximum  positive  strain.  Figure  6  shows  good  agreement  of  the 
radial  strains  at  170  mm  for  all  times.  Also,  it  appears  for  t  >  250  /xs  that  the  ADINA  curve  is 
beginning  to  lag  behind  the  EPIC-2  curve.  Figure  7  shows  good  agreement  of  the  radial  strains 
at  240  mm.  Figure  8  shows  good  agreement  of  the  normal  accelerations  at  106  mm  for  t  <  75 
/IS  and  poor  agreement  for  t  >  75  /is.  The  ADINA  positive  and  negative  peak  accelerations  are 
much  larger  than  the  EPIC-2  peak  values.  Figure  0  shows  good  agreement  of  the  normal 
accelerations  at  244  mm  for  t  <  150  /is  and  poor  for  t  >  150  /is.  Figures  10  and  11  show  the 
ADINA  and  EPIC-2  primary  shock  spectra  obtained  from  the  normal  acceleration  histories 
shown  in  Figures  8  and  9.  A  recursive  filter  technique^  is  used  to  obtain  these  spectra.  Figure  10 
shows  excellent  agreement  of  the  spectra  for  frequency  f  <  50  kHz  and  fair  agreement  for  f  >  50 
kHz.  Figure  11  shows  excellent  agreement  of  the  spectra  for  f  <  10  kHz,  good  agreement  for  10 
kHz  <  f  <  20  kHz,  and  poor  agreement  for  f  >  20  kHz. 


Figure  3.  EPIC-2  and  ADINA  normal  displacement  histories  at  106  mm  for  v 


P 


=  1508  m/s. 
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Figure  4.  EPIC-2  and  ADINA  normal  displacement  histories  at  244  mm  for  v 


P 


=  1508  m/s. 


Figure  5.  EPIC-2  and  ADINA  radial  strain  histories  at  100  mm  for  Vp  =  1508  m/s. 
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Figure  6.  EPIC-2  and  ADINA  radial  strain  histories  at  170  mm  for  Vp  =  1508  m/s. 


Figure  7.  EPIC-2  and  ADINA  radial  strain  histories  at  240  mm  for  Vp  =  1508  m/s. 
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Figure  8.  EPIC-2  and  ADINA  normal  acceleration  histories  at  106  mm  for  v 


1000  10000  100000 
Frequency  (Hz) 


Figure  11.  EPIC-2  and  ADINA  primary  shock  spectra  at  244  mm  for  v  =  1508  m/s. 
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IV.  SUMMARY  AND  CONCLUSIONS 


By  modifying  the  EPIC-2  code  to  calculate  the  time  rate-of-change  of  the  target’s  linear 
momentum,  the  code  can  calculate  a  loading  history  of  an  impacting  projectile  on  the  target. 
The  loading  history  of  a  20  mm  fragment  simulator  projectile  for  an  impact  velocity  of  1508  m/s 
on  a  914  mm  x  70  mm  thick  circular  rolled  homogeneous  armor  plate  has  been  calculated  by  the 
EPIC-2  code  and  has  been  used  in  an  ADINA  finite  element  analysis  of  the  plate.  Comparisons  of 
the  EPIC-2  and  ADINA  response  calculations  show  that  the  EPIC-2  calculated  loading  history 
can  be  used  as  the  loading  input  for  finite  element  analysis  of  ballistic  shock. 
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3  Air  Force  Armanent  Laboratory 
ATTN: 

AFATL/DLJW 
AFATL/DLORD 
AFATL/DLYV 
Eglin  AFB,  FL  32542-5000 

2  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  R-13 
Code  U-11 

Silver  Spring,  MD  20910 

2  Commander 

Naval  Weapons  Center 
ATTN:  Code  3181 
Code  3261 

China  Lake,  CA  93555 

1  Orlanda  Technology,  Inc. 

ATTN:  D.  Matuska 
P.O.  Box  855 
Shalmir,  FL  32579 

1  Southwest  Research  Institute 
Dept,  of  Mechanical  Sciences 
ATTN:  Dr.  C.  Anderson 
8500  Culebra  Road 
San  Antonio,  TX  78228 


No.  of 

Copies  Organisation 

1  AAI  Corporation 

ATTN:  R.  L.  Kachinski 
P.O.  Box  6767 
Baltimore,  MD  21204 

1  ADINA  R  &  D,  Inc. 

71  Elton  Avenue 
Watertown,  MA  02172 

1  General  Dynamics  Land  System  Division 
ATTN:  Dr.  P.  Kersten 
P.O.  Box  2074 
Warren,  MI  48090-2074 

1  Honeywell,  Inc. 

ATTN:  Dr.  G.  Johnson 
7225  Northland  Drive 
Brooklyn  Park,  MN  55428 

1  University  of  Denver 
Denver  Research  Institute 
ATTN:  R.  F.  Recht 
2390  S.  University  Blvd. 

Denver,  CO  80210 

1  Computational  Mechanics  Associates 
8600  LaSalle  Road,  Suite  614 
Baltimore,  MD  21204 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  laboratory  undartakaa  a  continuing  effort  to  Improva  tha  quality  of  the  reports  It  publishes.  Your 
eommants/ansHars  below  will  aid  us  In  our  efforts. 

1.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  Ir.terest  for  which 
the  report  will  be  used.)  _ 


2,  How,  specifically,  la  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas, 
etc.)  _ 


3.  Has  the  Information  in  this  report  lad  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate.  _ 


4.  Caneral  Comments.  What  do  you  thin)(  should  be  changed  to  improva  future  reports?  (Indicate  changes  to 
organization,  technical  content,  format,  etc.)  _ 


BRL  Report  Number  _  Division  Symbol  _ 

Check  here  if  desire  to  be  removed  from  distribution  list. 

Check  here  for  address  change.  _ 

Current  address:  Organization  _ 
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